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The retardation and temperature effects in two-body collisions are studied. The collision integral with 
retardation effects is obtained on the base of the Kadanoff- Baym equations for Green functions in 
a form with allowance for reaching the local equilibrium system. The collisional relaxation times of 
collective vibrations are calculated using both the transport approach and doorway state mechanism 
with hierarchy of particle-hole configurations in heated nuclei. The relaxation times of the kinetic 
(^ I method are rather slowly dependent on multipolarity of the Fermi surface distortion and mode of the 

CN ■ collective motion. The dependence of the relaxation times on temperature as well as on frequency of 

collective vibrations is considered and compared. It is shown that variations of the in-medium two- 
body cross-sections with energy lead to non-quadratic dependence of the collisional relaxation time 
00 ' both on temperature and on collective motion frequency. 

(N 

PACS: 05.20.H-w, 21.60.Ev, 21.65.H-f, 24.10.Pa, 24.30.Cz 
> 

(N 

^ : 1 Introduction 

O . The damping of the collective excitations as well as transport coefficients for viscosity and heat conduc- 

tivity are strongly governed by the particle collisions. The relaxation time method is widely used as the 
simplest and rather accurate approach for simulation of the collisional relaxation rate Ac oc 1/r, where r 
is the so-called relaxation time [jT], ^, |3|]. Relaxation time method can be applied to description of the decay 
rate of arbitrary mode of motion, but an explicit form of the relaxation time depends on specific features of 

y ■ the mode. In this contribution, the collisional relaxation times responsible for the width of the collective 

pj ! vibrations are studied. 

>: 

^ ■ 2 Semiclassical kinetic equation approach 

The collisional relaxation times can be calculated using the collision integral of the transport equation. 
In studies of the damping widths of collective excitation in the Fermi liquid, they are determined by the 
coefficients r^ of the multipole expansion of the total number J\f^^' (p) of the collisions in the direction 
P = p/p of the momentum space [^ Q] 

POO ,-A^ -^ 

ArW(p)^/ d6jW(p,e)= E E %Yimip). (1) 

Here, J^^'{p, e) are the linearized collision integrals 

jW (p, e) = J(±) (p, r, t) = ( Jp(p, r, t) ± J„(p, r, t))/2, (2) 

where the signs (+) and (— ) stand for isoscalar and isovector modes of vibrations, and the subscripts 
p and n stand for protons and neutrons, respectively; Jq,(p, r, t) is a collision integral in phase space 
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(p, r), when a nucleon of the sort a = (p or n) with momentum p is scattered; e is nucleon energy. The 
colhsion integrals are linearized with respect to the dynamical component of the phase space distribution 
function 6fa{p, r, t). The values £q determine the initial components of the multipole expansion of the 
total number of the collisions. The functions v^^^' are the partial components of the energy-integrated 
distribution function 5/(±)(p, r, t) = {6fp ± 5fn)/2 = <5/(±)(p, e, r, t), 

rdedf^^\p,e,r,t)=J2 E <p'£ir,t)YUp), 0) 

■'^ e>Om=-£ 

where Y^^ (p) is the spherical harmonic function. In approximation of a small difference in the chemical 
potentials for protons and neutrons and assuming fp = fn = f, where / = /(p, r) is the equilibrium 
distribution function, the dynamical distortions 5/*^^^ (p, r, t) of the phase space distribution functions are 
solutions of the linearized Landau-Vlasov equation 

-^^+ ^ ^ i ^ jf = J^^>ip,r,t), (4) 

at m or or op or op 

where 517^^'^ = 5U^^\p,r,t) is the Wigner transform of the variation of the self-consistent potential 
with respect to the equilibrium value e*^^-*. In the nuclear interior the mean field variation 6U^^' can be 
expressed in terms of the Landau interaction amplitude F^^^ (p, p/) as 

where Np = 2pF m* /{g tt"^ 'h?), pp is the Fermi momentum, m* is the effective mass of nucleon and g 
is the spin degeneracy factor. The quantity F^^^ (p, p/) is usually parameterized in terms of the Landau 
constants Fq and F^ as 

fW(p,p/) = Fo(±)+Fi(^)(p-p'). (6) 

In the isoscalar case, the Landau constants are related to the incompressibility modulus K [^ of matter 
and the effective mass m* [^] by 

K = 6/x(l + Fq^^^), m* = m (l + F^^+Vs) • (7) 

Here m is the mass of free nucleon and p is the chemical potential. We have that p ^ ep = p'\'/2m* 
for T <^ ep, where ep is the Fermi energy and T is the temperature. In the isovector case, the Landau 
parameter Fq' is related to the nuclear symmetry energy 6symm- Namely [0, ||], 

bsymm = ^fi{l+Fl^-^). (8) 

The quantities rj in Eq.(|l|) can be considered as the partial collective relaxation times because they 
determine a coUisional contribution to the damping widths resulting from the two-body collisions in the 
layer of the momentum space with multipolarity i, 

j^ dej dnpJ^^\p,e)Yeo{p)/ J^ de j d%5f^^'^Y,^{p) . (9) 
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These times are proportional to the relaxation times Tc defining the damping widths Tc (L) of the 
isoscalar and the isovector vibrations with frequency w [Q, ^ ^, |T^ in regime of rare collisions with 



iJTc ^ 1 in the Fermi liquid. In particular, the coUisional damping widths of giant resonances with 
dipole (L = 1) and quadrupole (L = 2) multipolarities resemble the widths in the relaxation rate approach 

rW(L) = n/rW(L), r(-)(L = l)=rti, ri+)(L = 2)=r£), (10) 



in the case when nuclear fluid dynamical model with relaxation is used ||4j, ^, |10|, |11|, [12p. The collisional 
damping width[^] of zero sound in the Fermi liquid with its relative velocity 5*^ ~ 1 is also given by 
Eq.(|l^) but with the use of the t^^'^ oc T2 for Tc (L). The time tJ;^2 at a; = is the thermal relaxation 
time determining the viscosity coefficient of the Fermi liquid [13]. 

The variations of the mean field and of the dynamical component of the phase-space distribution func- 
tion change rapidly in the systems with high frequency collective vibrations. This leads to the memory 
(retardation, i.e. non-Markovian) effects in the collision term. There are different expressions for memory- 



dependent collision integral in the Fermi liquid ([14]- [20]). 



The non-Markovian collision term of the semiclassical Landau- Vlasov equation was studied in Bom 
approximation with the use of the Kadanoff- Baym equations for the Green functions in Refs. [ pTJ [T9| ]. 
The one-component Fermi liquid was considered with the periodic time variation of the nonequilibrium 
distribution function df = 6fn = Sfp, 6f oc exp(— iwt). As a result, the linearized collision integral 
consists of two components (see Eqs.(42),(43) and (45),(46) of Ref.[|T9[]), i.e.. 



J(p,r,t)= J«(p,r,t) + J(2)(p,r,i), (11) 

where the components j(^)(p,r,t) and j(^)(p,r,f) are determined by the variations of the distribution 
function and the mean field, respectively, and 

J^^Hp,r,t) =2 J ^^^^0^Wi{p.mAp)B(^Hp,r,t). (12) 

Here, l^({pj}) = {da/di^)4:{27rh)'^/m'^ is the probability of two-body collisions with the initial momenta 
Pi = P)P2 and final ones P3,P4, {i = 1^4); da/dQ is in-medium differential cross-section (in Born 
approximation); 

S«(p,r,t) = ^ 6h{t) ^MM. [<5+(Ae- + no;) +5_(Ae-- 7i^)], 
k=i '^ Jk 

(13) 

i?(2)(p,r,t) = g({/^.})^M!)){[5^(Ae-+ n^) - 5+(A6-)] - [5_(A6-- ^o;) - <5_(A6-)]}, 

where Jk = /(pfci r); dQ{{fj})/dfk are the derivatives of the is the Pauli blocking factor Q with respect 
to/fc, 

QiUj}) = (1 - /l)(l - /2)/3/4 - /l/2(l - /3)(1 - h). (14) 

The ej = e(pj,r) is the equilibrium single-particle energy and dUj the variation of the mean field, and 
Ae = ei + 62 — es — €4, A{6U) = 5Ui + 6U2 — SUs — 5U4, Ap = pi + P2 — Ps — P4. The equilibrium 
distribution function /^ = f{pk , r) depends on the equilibrium single-particle energy e^ = e{pk , r) : 
fk = fi^k)- It equals the Fermi function evaluated at the temperature T, /(efc) = l/[l+exp((eA; — fi)/T)]. 
The nonequilibrium component 6f of the distribution function can be presented in the the form 

6f{p,r,t)= -u{p,T,t)^. (15) 



With the use of this relation the expression for quantity B^^' can be transformed to the form 
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k=l 



Ai^Qiifj})^ [5+{A€ + huj) + 5^{Ae-hu)] - 5B^^\ (16) 

anuj 



where Av = vi + 1^2 — V3 — 1^4, Vk = ^(Pfe, r, t) and 
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E Q({/,})['^+(Ae- + ncu) + (5_(Ae--7iw)]|^. (17) 



The first component in the Eq.(|17D determines a probability flux of colliding particles which is con- 
nected with possibility of variation of the energy ek when the values of other energies (e^yfc and huj) are 
fixed. This term should be equal zero because of fixing the total energy in initial or final states and therefore 
it does not contribute to the total number of the collisions J\f^^^ , Eq.(|l|). The last statement can be easily 
verified by direct calculation of this contribution to the J\f^^' with the use of the procedure proposed by 
Abrikosov and Khalatnikov (see Eqs.(25)-(|30|)) for evaluation of the manifold energy integrals. A relative 
dynamical component Uk of the distribution function is slowly dependent on energy and it can be consid- 
ered (at least for low temperatures T ^ e^) as a function of the momentum direction rather than of the 
momentum: u^ = i^(pA;,r,t) = i^(pfc, ei?, r, t). Therefore the second component in the Eq.(|l7|) is also 
negligible and the term 6B^-^^ in the Eq.(p^) should be rejected, 

(5-B(^) = 0. (18) 

Note that the generalized functions 6j^, 6^ appearing in Eqs. ([l3|), ([iq) and (|l7|) include also integral 
contribution, 

6^{x) = ^ r dr e— = f -^ = hix) - -^T'l-), 6.{x) = dl{x), (19) 

^vr J -00 2tt X + lU 2 2m x 

where 6{x) is the delta function and the symbol V denotes the principal value of integral contribution. 
The integral terms of the 6±, corresponding to virtual transitions, are usually omitted in the J because 



they assumed to be included by renormalizing the interactions between particles [21 1. This corresponds 



to substituting 5{x)/2 for 5± in Eqs. (|13|), i.e., to taking into account real transitions with conservation 
of energy. The shift in energy Ae by 7i a; in the arguments of the (5-functions of the expressions for the 
collision integral agrees with the interpretation of the collisions in the presence of the collective excitations 



proposed by Landau [22]. According to this interpretation, high-frequency oscillations in Fermi liquid can 
be considered as phonons, that are absorbed and created at the two-particle collisions. 

In the one-component Fermi liquid the nonequilibrium distribution function (5/(p, r, t) = /(p, r, t) — 
/(e(p, r)) is a solution of the linearized Landau- Vlasov equation in the form 

^ + iL(p «„7=J, (20) 

ot m or 

Here, 6f is a linear deviation of the distribution function from its local equilibrium value fi,e., where a 
function fi^. is equal to the the Fermi function /(e) evaluated with actual one-particle energy e = e + 6U, 

/,e. = /(e(p,r,t)), 

Sf = df-^6U = f{p,r,t)-fu,= -x%, X = ^ + SU, (21) 

ae ae 



where relationship for 5U has the form of the Eq.(|5]) but with the interactions and distribution function for 
one-component Fermi liquid. 



According to Eq.(20), a Fermi system tends to the local equilibrium (when d6f/dt = 0) if the collision 
integral is a functional of the 6f, J = ^{6f). The collision integral given by the expressions ([rT|)-([l3|) 
and (16), (18) has the following general form: 5 J = 6J^^'{6f) + 5J^^'{SU). Therefore it does not lead 



to local equilibrium of a system. The condition of the existence of the local equilibrium of a system is 



general property of the Landau- Vlasov equation in the Fermi liquid [23, E4] at d6f/dt = 0. Therefore 



the Born approximation (|11|)-(|13|) for collision integral is poor approach without additional modification, 
and a revision of the derivation method of the collision integral expression is needed. 

It should be initially noted that the foregoing relationships for collision integral are obtained with the 
use of the perturbation theory in nearly nonhomogeneous systems with week interaction between particles. 
A week interaction can not change rapidly the trajectory of the particle and due to this it can not lead 
to rapid variations of the distribution function. It means that retardation effects are overestimated in the 
expression for collision integral in Bom approximation where it was assumed that distribution function was 
varied very quickly during all possible interval of the time changing (— oo < t' < t). Consequently, the 
collision integral given by Eqs. ([TT|)-([T3|), ([l6|), (p^) can be in fact correct in the case of small retardation, 
i.e., for small values of the fku. 

With this in mind, we replace the derivatives of the form dd^{Ae + huj)/dhio and d6-{Ae — 
fiu))/dTiujm the Eq.(|l|) by the finite differences ((^+(Ae + ?iw) - 5+{Ae))/nu) wA ((5_(Ae) - (5_(Ae - 
Uu}))/fiuj, respectively. Then we combine the resulting expression together with contribution B'^'^' arising 
from mean-field variation and obtain the linearized collision integral for one-component Fermi liquid in 
the following form 

J{p,r,t)=J (2^;,)6 ^({p.})'^(Ap)AxQ . (22) 

With the use of the algebraic relation [^3]] 

[(1 - /i)(l - /2)/3/4 - /i/2(l - /3)(1 - h) exp (^^j]5(Ae ± hu;) = 0, (23) 



the Eq. ( [22| ) can be presented as 



I ^^0^W{{p.})5iAp)Axfif2{l - /3)(1 - h) mnu;,T) - <l>i-nu;,T)] , (24) 



where <^{huj,T) = 6{Ae + huj)[exp{-huj/T) - l]/nuj. 



The collision integral of the form (g2|) or (g4|) provides driving distribution function towards its local 
equilibrium value because it depends on the variation df, J = J{6f). This behaviour is in line with general 
properties of the Landau- Vlasov equation in the Fermi-liquid [ |23| , |24] ] at d6f/dt = 0. The expressions 



(22), (24) depend only on the occupation probability 7-'2p2/i = /i/2(l — /3)(1 ~ Ja) of the 2p-2h states 
in the phase space. This fact leads to interpretation of the collisional damping with Unearized collision 
term as the relaxation process due to the coupling of one-particle and one-hole states to more complicated 
2p — 2h configurations. 



The form of the collision integral (24) in the Markovian limit (lu —>■ 0) coincides with the standard 



expression for the collision integral in Fermi-liquid without retardation effects [ ]23| , |24[ ] because in this case 
the term in square brackets of Eq.(|^) tends to the value —6{Ae)/T. 



The equation (22) for some special explicit form of the quantity Xj was used at first in Refs. 1 16, 25, 26]. 



The derivation of the collision integral ( |22D is performed in Ref. | |14| ] within framework of the extended 
time-dependent Hartree-Fock model. The expressions for the distortion functions Xj corresponding to a 
perturbation approach on collision term and including the amplitudes of the random phase approximation 
were used in this method. 



The expression for the collision integral in two-component Fermi-system is obtained from Eq.(g4|) in 
the same manner as done in Ref.[Q] under the assumption that chemical potentials and the equilibrium 
distribution functions are the same for protons and neutrons. 

The analytical expressions for partial collective relaxation times r^ ( Eq.(^ ) can be obtained in low- 
temperature and low-frequency limits (T, huj <C ep)- In this case the momentum integrals are calculated 
using the Abrikosov- Khalatnikov procedure[23, |^, 28, 29] which is based on the assumption that par- 
ticles are scattered near Fermi surface with the momentum values pi approximately equal to the Fermi 
momentum pp. In this case the probability TV ({pj}) of two-body collisions can be taken as a function of 
two scattering angles (p and 6, where (p is the angle between the momenta pi and p2, and 9 is the angle 
between the (P1P2) and (P3P4) planes; that is 



COSCp = (,P2-Plj, 

COS e = [pi X P2] • [P3 X P4] 



iPl XP2JIIIP3 Xp4j|. 



(25) 



It allows to separate the angular and the energy integrations in the collision integral at arbitrary scattering 



angle[29]. The integrals with respect to momenta in expression for the collision integral are calculated 



employing the transformation [23, 27, 29] 



(ip2dp3iip4<5(Ap)(. 



*3 

m 



dO 



cos 



A I d,f / de2de2.dei{...). 

cos t; Jo jo jo jo 



(26) 



Here, ip is the azimuthal angle of the momentum p2 in the coordinate system with the Z axes along pi. 



The integration with respect to the azimuthal angle ip is performed by the relation[|28|] 

/ TT^lrniPj) =YimiPl)PliPjPl), 

Jo ^vr 
where Pi is a Legendre polynomial, and 

(PsPi) = cos2(0/2)+sin2((/)/2)cose, 
(p4Pi) = cos2((/)/2) -sin2(0/2)cos6'. 



(27) 



(28) 



00 

+00 



To perform over energies in the collision integral the following expressions are used [ 23 ] 

/■+00 
^lyiy) = I dxidx2...dxyn{xi)n{x2) ...n{xy)5{xi+ X2 + .■■+Xu) = 

/+00 
dtn{t - y)Iu-i{t), 
-00 

where n{x) = 1/(1 + exp (x)), n{x) + n{—x) = 1 and 

1 y^ + n^ 



(29) 



hiy) 



2 1 + exp {-y) 



, h{y) 



1 y{y' + A7r') 
6(l-exp(-y))' 



(30) 



Finally we get the following relation for relaxation times using the collisional integral of the form given 
by Eqj 



n 



.(±) 



TZ{uj,T) 



< aL<^y' > + < a'^) ' > 



(±) 



pn - 



(31) 



where a'^^ = ((T^„ + u' )/2; a'--, = dajj'/dQ is in-medium differential cross-section for scattering of the 
nucleons j and / ( here, j = n or p, and similarly j' = p or n). The quantity 7l{uj, T) result from the 
energy integrations and has the following form in low-temperature and low-frequency limits (T, huo <^ ep) 
in the approximation m* ~ m 

2 m 

The symbol < ... > in Eq.(|3ll) denotes averaging over angles of the relative momenta of the colliding 
particles, 



7^(w,^) 



,_^{{hcof + {2nTf}. 



(32) 



vr Jo 



#sin(0/2) / de{ 



(33) 



The functions <1>^ in (^Tj) define the angular constraints on nucleon scattering within the distorted 
layers of the Fermi surface with multipolarity I: 



$ 



(±) 



1 ± P,(cos 4') - PeiiPsPi)) T PiiiPAPi)), 



(34) 
where the scalar products (psPi) and (p4Pi) are given by Eq.(|28|). It follows ^\Vq{(J), 9) = ^\Ji{4>, 9) = 



$ 



i=oy 



,(+), 



0. These relations lead to possibility of the two-body collisions in layers of the Fermi surface 



(+) 



^(-) 



1 for the 



distortion with multipolarity beginning with the value £q = 2 in the isoscalar case and 

isovector vibrations. As a result, the isovector dipole relaxation time t^^{ has a finite value, that means a 



nonconservation of the isovector current in the presence of n — p collisions [|30|]. 

Due to the momentum conservation and conditions pi= pp, the angle 9 agrees with the scattering angle 
in the center-of-mass reference frame of two nucleons. The angle 4> defines the magnitudes of the relative 
momenta kj = (p2 — pi)/2 and kj = (p4 — p3)/2 before and after collision, respectively. The value of 
total momentum, P = pi + p2, also depends on a magnitude of the (j). We have 



kjkj = k cos ( 



k' 



k 



k) = pI W (0/2) , P^ = Apj, cos^ (0/2) . 



(35) 



Therefore the relative kinetic energy Erei of two nucleons as well as the energy of centrum mass motion 
Ecm are dependent on angle 



Erei = k^/m = 2ep sin2(0/2), Ecm = P"^ /2m = 2ep cos^{(j)/2) 



(36) 



and the total energy Etot = Erei + Ecm holds only fixed, Etot = 2ep. Therefore the in-medium differential 
cross-sections a', „ of two nucleon scattering depend on the relative momenta k^ and kj at fixed total 
energy rather then at fixed relative kinetic energy Erei, because the magnitude of Erei changes with angle 
between colliding particles. The transfer momenta q = kj — kj = p3 — pi and q' = — (kj + k/-) = 
Pi — P4 for scattering due to direct and exchange interactions respectively are also functions of the and 
9:q = 2k{(j)) sm{9/2) and q' = 2A;(0) cos(6l/2). 

Now we estimate the collisional relaxation times in the case of the isotropic scattering with independent 
of energy the angle-integrated cross sections ajji. Performing angular integration in (pl|) with the use of 



Eqs.(^ and (g^) we find that 1/r 



(±) 
i<e. 



(±) 



Oand 



n 



1 



.(±) 



a 



(±) 



{hLj/2lTf + T^ 



1 



a 



8m 



(±) 3h' 






np 



(37) 



Q 



1 



-1 



d 



{-) 



l-(-l)^ 



d' 



(+) 



d 



(-) 



Q=o = Q=i = 0, 



2^ + 1 ' '^ 2£+l 

where aav = [o"pp + <7nn + 2(Tnp] /4 is the in-medium spin-isospin averaged nucleon-nucleon cross section. 
The magnitude of the in-medium cross section ajji is taken usually proportional to the value of the free 

jo^^i, so that the parameter a^ can be rewritten in 



space cross section a--, with a factor F 
the form 



^n 



n 



a 



(±) 



a 



(±) 



IF, 



a 



(±) 



4.18/ 



Q + l.Sd 



(±) 



MeV. 



(38) 



(free) 



5 /?7z^ are adopted|[T^, |lO|]; they correspond to the 



Here, the values aJv'^'^ = 3.75 fm? and anp 
free space cross sections near Fermi energy. 

The relative relaxation times rj /'^1=2 ^i*^^ '^^^ f''^^ space cross sections are shown on Fig. 1 in re- 
lation to the multipolarity i of the distorted layers of the Fermi surface which are accessible to particle 
collisions. Solid and dashed lines connect the values which correspond to isoscalar and isovector modes 
of vibrations respectively. The magnitudes of the relaxation times are different for isoscalar and isovector 
modes of vibrations and they are dependent on the multipolarity I. As seen from the Fig.l, the coUisional 
relaxation times rather slowly vary with multipolarity £ and with collective motion mode at isotropic scat- 
tering with energy independent free cross sections. In particular, parameters ai , which define relaxation 



times by the Eq.(|37|), take the same value at 



a 






1.1; a/a 



(+) 

1=2 



0.8, 5£) 



/a 



(+) 

(.=2 



~(±) 



= a 



4.18 MeV, and 5/5^ ~ 0.9 
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Fig. 1 : The relative relaxation times r^ /rg ver- 
sus multipolarity £.Solid and dashed lines con- 
nect the values which correspond to isoscalar and 
isovector modes of vibrations respectively. 

The dependence of the viscosity relaxation time r)^2 (^ 



T,MeV 

Fig.2: The relaxation time T2 in dependence of 
T athio = 0m units of 10"^^ sec. 



0, T) on the temperature is shown in Fig.2. 
The value of 0)^2 i^ used from Eq.([38[). The temperature dependence arises from smearing out the equilib- 



,.(+) 



rium distribution function near the Fermi momentum in heated nuclei.The coUisional rates l/r) ' given by 



(i) 



Eqs.(pll) and (37) are quadratic both in temperature and in frequency with the same relationship between 
the components much as in the zero sound attenuation factor of heated Fermi liquid within the Landau 



(±) 



prescription [22, H, 16, 14 1. The relaxation times r^ depend on frequency uj due to the memory effects 



in the collision integral. 



3 Doorway state mechanism in heated nuclei 

The relationship ( [Iq| ) gives the possibility to evaluate the relaxation time in system with weak damping in 
an independent way from decay rates Ac = rj; /fi. We adopt the following physical notion: Ac = 
Tc /fi is the spreading decay rate of the initial state \v\ ) to the final state \vr ') within first- order 
approximation of the time-dependent perturbation theory, as given by the golden rule 

AW = l/rcW = ^^if(^pf), (39) 



n 



where p V is the density of the available final states. 



The quantity |M(='=) p is the mean square matrix element for transitions due to residual interaction Vn 



\Mi±)\'^ = \{uf\Vres\y^^)\\ (40) 

where the line over symbols denotes an average over final states ^ [^. 

The initial state should describe giant collective vibration in heated system at given temperature T. It 
is taken as a mixture of a collective state {GR) and a thermal state which is approximated by uncorrelated 
superposition states of mp — mh configuration with m excited particles and holes corresponding to the 
most probable number of excitons n = 2m at given temperature T. The excitation energy of the system 
is the sum of collective energy Ulo and thermal excitation energy U = ne with e = 7r^T/(12 In 2) for the 



average excitation energy per thermal exciton [ |34| ] : 

\v^) = \{GR], {mpmh}), E = huj + U, U = aT"^, n = 2m = 2gTln2, (41) 



where the expression for n is taken from [Q]; a = TT^g/6. The quantity g is the single nucleon state 
density at the Fermi surface and the same values of g are taken for neutrons and protons. 

Next we accept common feature that giant resonance state [GR) is formed by coherent superposition of 
many (predominantly correlated) one-particle one-hole configurations and due to this fact wave function 
of initial state can be presented as the sum of wave functions \{{m + l)p{m + l)h}) = |{nj}; ki) of 
incoherent (tti + l)p — (tti + l)/i configurations with n^ = 2 + n excitons, fc, stands for other quantum 
numbers 

Wi) = T.Ck:\{^^}■,ki), (42) 

where the quantity C^' defines the magnitude of the admixture of different components of quasiparticle 
eigenstates. 

Because of two-body character of the residual interaction Vres, the final state can consist of configu- 
rations with nj = Hi, Hi ±2 excitons. The averaged squared matrix elements |Mp of the transitions to 
states with fixed number of excitons can be rewritten as 






where M.'^{ni —>■ nf,E) = |({n/}|V^es|{?^i})P i^ effective mean square matrix element for transition 
between incoherent particle-hole states. 

"To simplify the presentation, we will omit in the following the superscript (±) and include them only when it is necessary to 
avoid confusion. 



This transformation is performed by the use of the following assumption and properties: 
i) The compensation of the binary products of the matrix elements coupling together incoherent exciton 
states with different values ki is assumed to take place due to very complicated character of the final state; 
ii) Approximate normalization of the factors C^, is used, J2k- l^'^'P — 1; 

iii) The mean square matrix elements for transitions between incoherent exciton states with different values 
of the numbers ki, kf are taken as equal to the same magnitude 7W^(ni — > Uf, E) which is dependent only 
on numbers of excitons and the total excitation energy. We also assume that effective mean square matrix 
elements Jvp{ni — > nj, E) for interactions between different kinds of nucleons are equal in magnitude [|^, 



33,34,35]. 



With the use of (|43|), the coUisional relaxation rate Ac (Eq.(^) coincides with the particle interactions 
rate of the exciton model starting from the ni configuration] 34]. The relation ( p9| ) for the collisional 
relaxation time Tc = Tc{uj, T) is 



— = 2TTM^{n, 



E = hLO + U, U = aT, n 



n + 2,E) pciE) + 2^X2(n< = n, E) pa{E), 



(44) 



hT, 



'k'^q/Q, 6 = 0.843a, 



when processes of creation and of annihilation of the particle-hole pairs are included. The matrix elements 
for both processes are taken to be determined by the number excitons n< in the simplest state [^ and pc 
(Pa) is the density of the final accessible states corresponding to the pair creation (annihilation). 

The transitions to final configuration with nj=nj + 2 = n + 4 dominate at low excitation energies. 



Using the simplest expression within the exciton model [34] for density of final accessible states, Pc{E) 
{g^/2){E'^/{ni + 1)), the Eq.® is given by 

{hoj + aT'^f 



n 



Trg^M' 



.n< 



n + 2,E = huj + U)- 



(45) 



Tc{uj,T) ■■" " " '■^■^ '^ ' ' ■'" ' " ' 3 + bT 

According to the exciton model studies ([33]-[40] the effective mean square matrix elements Ai'^{ni,E) 
is energy-independent at low excitation energies and it is inversely proportional to energy at higher excita- 
tions. The energy-independent estimation AA."^ was obtained with the use of the Fermi gas model as[|^, |^] 
M."^ = Km/-^^, K,m — 15.3 MeV^, where A is the mass number. The behaviour of collisional relaxation, 
as given by Eq.(p3|), with such magnitude of the mean square matrix element agrees with estimation (|^) 
based on kinetic equation approach at low temperatures T <C fiuj. 

There are different estimates for the mean square matrix element with dependence on energy and 
number of excitons [33, 39, BO]. The fulfillment of the condition of equiprobability of all particle-hole 



configurations is assumed in most of them and therefore they can not be used in the considered case of 
collective (predominantly Ip-lh) state overlapped with temperature-fixed background particle-hole states. 
The expression for M.'^{n, E) without assumption on a uniform sharing of the excitation energy E into n 
excitons was proposed in Ref . [ |40| ] : 

n< + 1 Kb 



M' 



n. 



:,E) 



(46) 



4 A^E' 

where quantity Kb is not changed with E and n but can be dependent on numbers of protons and neutrons, 
Kb = 190 MeV^. If this value of Jvi"^ is employed as the squared intronuclear matrix element, the 
collisional relaxation time Tc is a linear function of the collective energy tiuj and thermal energy U. The 
corresponding expressions (Eqs.(p5|) and (46)) for the relaxation time have the same form as that one 
obtained within test particle approach, when collisions were simulated by modeling s-wave scattering 

^: 



between pseudoparticles [|41 



n 



1 



Tc{u},T) ae 



{Tioj + U) , 



/Csvr 



Or 



{g/Af 



(47) 
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The relaxation times given by Eqs.(^7|) and ( p7| ) can agree together at zero temperature if the magnitude 
of Kb is equal to the value /Co = {huj/a)[A/{gT:))'^ . Here, /Co = 70.9/iw/a ~ 220 MeV^ for giant 
isovector dipole resonances in heavy nuclei, when huj ~ 13 MeV, g = A/13 and a = a = 4.18 MeV. 
This value of Kb is rather close to the Kb = 190 MeV^. It means that in cold nuclei the relaxation times 
for the GDR within doorway state mechanism are not too different from those obtained within the transport 
approach. 

The dependence of the coUisional relaxation times Tc{uj, T) on temperature and energy fiuj is demon- 
strated on Figs. 3-6. Solid and dashed lines correspond to the relaxation times Tc{uj, T) within doorway 
state mechanism with the mean square matrix elements A4^(n<,£') oc 1/E and M^{n^,E) = const 
respectively. Dot-dashed lines correspond to the collisional relaxation times Tpl within the framework of 



„(-) 



the transport approach with the value a\ determined by free n — p cross-section ( see, Eq.(|38|)). The 
factors Km and Kb of the mean square matrix elements are fixed from the condition of the coincidence 



of the relaxation times Tc{u! = Egdr/^, T = 0) and T^ 



(-), 



LJ 



Egdr/^,T = 0) in cold nuclei at a 



frequency corresponding to the GDR energy. The magnitude of this energy is taken as equals to the GDR 
energy in "^^^Pb: Egdr = 13.43 MeV. The values of the relaxation times are given in units of 10~^^ sec. 
Figures 3, 4 show relaxation times at w = (Fig.3) and u = Egdr/^ (Fig-4) in relation to the tem- 
perature. 



o 

CO 





T, MeV 

Fig.3: The relaxation time r in dependence of T 
athuj = in units of 10^^^ sec. Solid and dashed 
lines denote r within dooway state mechanism with 
A/P oc 1/E and AP oc const. Dot-dashed line is r 
within the framework of the transport approach. 

Figures 5, 6 demonstrate dependence of the relaxation times on energy e 
and heated (Fig.4, T = 2 MeV) nuclei. 



12 3 4 5 

T, MeV 

Fig.4: The relaxation time in dependence of T at 
huj = Egdr in units of 10^^^ sec. Notations are 
the same as in Fig. 3. 



hoj in cold (Fig.3, T = 0) 
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8, MeV 



10 15 

£,MeV 



r 
20 



Fig.5: The relaxation time in dependence of e = Fig.6: The relaxation time in dependence of e = 
?ia; at T = in units of 10^^^ sec. Notations are tku atT = 2 in units of 10^^^ sec. Notations are 
the same as in Fig. 3. the same as in Fig. 3. 



The collective relaxation times in heated nuclei can be approximately presented by the expression 



l/r,{iO,T) = qi[J' + q2T']y[q3 + q^T]\ 



217, 



(48) 



where qj are some constants and the exponents are the functions of frequency and temperature and they 
are /3 = 2, 7 = 1 , (5 = in the transport method; the /3, 7 are changed from 2 to 1 and 6 varies from 1 to 
with growing of the excitation energy in the doorway state approach with allowance for pair creation. 



4 Results and conclusions 

The retardation and temperature effects in two-body collisions in heated Fermi-systems were studied. An 
expression for non-Markovian collision integral of the Landau- Vlasov transport equation was obtained in 
a form which allows for reaching the local equilibrium in system. It was found in a small retardation limit 
on the base of the Kadanoff- Baym equations for Green functions. 

The expressions for collisional relaxation times of the collective vibration in heated nuclei are derived 
with the use of the non-Markovian collision integral as well as of the decay rates of exciton model. The 
relaxation times depend on frequency of the collective vibrations and the temperature. The temperature 
dependence arises from smearing out the equilibrium distribution function neai" the Fermi momentum in 
heated nuclei. The frequency dependence results from the retardation (memory) effects in the collisions. 
Analytical expressions for relaxation times of the isoscalar and isovector modes of the collective motion 
are derived in the case of the energy independent isotropic cross-sections in the two-body collisions. The 
relaxation times rather slowly vary with multipolarity of the Fermi surface distortions governed by collec- 
tive motion and two-body collisions. It gives possibility to use approximately the relaxation time ansatz 
for collision integral. The relaxation times depends on type of collective motion mode hke the lifetimes of 



the particle-hole configurations in two-component exciton model of the Ref.[]32|]. 

New approach for calculation of the collision relaxation time in heated nuclei are proposed using the 
formulae for the transition rates of the particle-particle transition between thermal state with collective 
vibrations and incoherent particle-hole configurations. This method leads to the same results as the trans- 
port approach in the case of low temperatures and energy independent mean square matrix element of 
interparticle collisions. It makes possible to take into account the energy dependence of the in-medium 
cross-sections in a simple phenomenological way by the use of the parametrization of the mean square 
matrix element 7W^ of interparticle collisions from exciton model of nuclear reactions. The dependence of 
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the matrix element M'^ ( i.e., the in-medium cross-section) on energy leads to non-quadratic dependence 
of the relaxation times on temperature and collective vibration frequency. 
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